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This work is aimed at demonstrating the possibility to construct new exactly-solvable stochastic 
systems by use of the extended supersymmetric quantum mechanics (N — 4 SUSY QM) formalism. 
A feature of the proposed approach consists in N = 4 SUSY QM the fact that probability densities 
and so obtained new potentials, which enter the Langevin equation, have a parametric freedom. The 
latter allows one to change the potentials form without changing the temporal behavior of the probability 
density. 



1 N = 4 SUSY QM 

Let us get started with a brief discussion of the N = 4 SUSY QM structure and con- 
structing within its framework isospectral Hamiltonians [H |2] . The N = A SUSY QM 



Hamiltonian has the following form: 
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H n „ = - 2 ^ + V ?(*) + 4 ] V 2 \x)) = -(p 2 + V^x) + afWx)) (1) 
where we have introduced (h — m — 1): 

a: KM-i"M + 5^g P) 

W(x) is a superpotential, are matrices which commute to each other and have the 
eigenvalues ±1 

o-W = <7 3 <g> 1, af ] = 1 ® a 3 (3) 
Supercharges Qi of the extended supersymmetric mechanics form the algebra: 

{Q h Q k } = 25 lk H; {Q u Q k } = 0;i,k = 1,2 (4) 

And admit the form: 

Q t = a ( l\p + iV^ +1 \x)), Qi = a ( ;\p - iV^ i+l \x)) (5) 

Constructing isospectral Hamiltonians within the iV = 4 SU SY QM is based on the 
fact that four Hamiltonians cast into one supermultiplet. Let us take as an initial one of 
the Hamiltonians 

#<ri,ff 2 = \{p - iviV a2 (x))(p + ioiV a2 (x)) + e = ~(p - ia 2 V ai (x))(p + ia 2 V ai (x)) +e 

(6) 
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where 

V a2 (x) = W'(x) + ±<r 2 W'\x)/W\x), (7) 

and e is the so-called factorization energy. In what follows the energy level will be counted 
from e. 

Consider an auxiliary equation: 

H%<p(z) = e<p(x) (8) 

Its general solution (p(x,e,c) is the linear combination of two independent solutions 
cp^(x, e) and cp^(x,e). The expression for W(x) in terms of (p(x, e, c) has the from 
of: 

W(x) = yln \l + \jdx' [cp{x', e, c)] 2<J1<J2 J (9) 

where A, two new parameters. 

For defmiteness we will consider <Ji = 02 = — 1. Then, relations in the above can be 
presented in the form of: 

d 2 1 W \<p(x,e,c)ibz(x,E)} 

dx 2 + ^ 2 {E - e) <p(x,e,c) 

H+ = Hz-£ 2 ln(l + XJdx' [tp(x', e, c)f) 
i/)$(x,E) = ^z(x,E) P^£A jdx' 2 <f(x',e,c)^z(x',E) ( n ) 

l+A J dx'[tp(x',e,c)] 2 x i 

x i 

In the case when (p(x,e,c) is a ground state eigenfunction of the initial Hamiltonian, 
wave function of ground state has a form: 

E ) = iV (12) 

l + A J dx' [<p(x',e)} 2 

Xi 

2 Construction of stochastic models associated with 

N = 4 SUSY QM. 

The Fokker-Planck equation is equivalent to the Langevin equation, however the Fokker- 
Planck is used more widely in physics, since it is formulated in more common for the 
probability density mf(x, Xq) language. The Fokker-Planck equation takes the form [31 H]: 



^- t mt(x,x ) = %-§^mt(x,x ) =F 4-^{x)mf(x, x 

+ x (13) 

m t=0 (x,x ) = (8(x-x )}] U±(x) = =F / dz$(z) 



with U±(x) to be the potential entering the Langevin equation. The Fokker-Planck equa- 
tion describes the stochastic dynamics of particles in potentials U + and LL = —U+. 
Substituting 

m?(x,x Q ) = exp{- [U ± {x) - U±(x )} /D}K ± (x,t) (14) 
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the Fokker-Planck equation transforms into the Schrodinger equation with imagine time: 



-D±K ± {x,t) = & + | [&{x] ±D®{x)]\K ± {x,t) = H ± K ± (x,t) 
K±(x,t) = (x \exp{-tH±/D}\x ) 



(15) 



in which the diffusion constant can be treated as the "Planck" constant. As in the previous 
case, let us consider H~ to the initial Hamiltonian. 



Vt\x) 



p 2 +(yt\x)) 2 -Dvt ) ' 



x) 



W'ix) - f ±hi\W'{x)\ = -D±\n<p(x,e,c) 
The force entering the corresponding Langevin equation: 

F(x) = 



(16) 



= D± In <p(x, e, c), UZ (x) = -D In <p(x, e, c) 



(17) 
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Then, in the same to the N = 2 SUSY QM way, there is the relation U + 
directly follows from the fact that Hz is obtained from Hz with changing V 2 '(x) to 
-Vt\x). 

Hence, there is the relation between the stochastic dynamics in a potential and that 
of in an inverse potential. 

Further consideration is based on an important property of N = 4 SUSY QM, which 
consists in having the symmetry of HZ\ under 0\ <-> cr 2 - It leads to: 

H-(x,p) = = iQ^Q^ = H+(x,p) (18) 

If the first equation points to the expression of H^ in terms of Qi(Qi), the second one 
is H^ in terms of supercharges Q2(Q2) ■ Though the supercharges of H^ and if* - * 1 are 
essentially different: 



Hz(x,p) = lQ?>Q ( 2 



p 2 +(v} +) (x)) 



DV{ +), (x) 



Vl + \x) 



£(W + %ln\W'{x)\)=D±]n 



ip(x,e,c) 



1+A J dx'[<fi(x',e,c)] 2 



dUZ 
dx 



(19) 



From equations ( fl9l) follows that at the quantum level the Hamiltonians Ht and Hz\ 
possesses the same spectrum and wave functions, thus the corresponding to them K(x, t) 
are the same. At the same time the corresponding to these Hamiltonians stochastic models 
are described by essentially different potentials, so their m t (x,Xo) are different. Further, 
the potential U^(x) has a nontrivial parametric dependence on A, that corresponds to 
having the family of stochastic models, the probability densities of which have the same 
temporal dependence. Having the parametric freedom allows one to change the potential 
form that is unexpected in the case. In particular, as it has been noticed in literature, 
physical quantities such as the time of passing a peak of potentials depend essentially on 
local changes of the potential barrier. Looking at Hf, we note that it follows from 



with changing Vi + \x) to —Vi^'(x), as it also takes place in iV = 2 supersymmetry. In 
its turn the potential of this stochastic system is U+(x) = —U^{x). The wave functions 
which are used to calculate K+(x,t) have the form of (fTTj) in the case. 

Getting back to the parametric dependence of mi t (x,x ) and m^ t (a;,xo), it should 
be noted that the normalization condition could lead to fixing the A. The same does 



'(+) 
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not happen in the case of m+ tt (x,xo), if instead of (p(x,e,c) we will use the normed 
wave function (say, the ground state wave function) of the initial Hamiltonian HZ- From 
definition: 



ip(x,£,c) 



1+A J dx'[<fi(x',e,c)] 2 



ip(x,e,c) 

in view of eq. ( TT2l) it follows that 



l+A J dx'[ip(x',£,c)f 

x ^-7—, E e-^+(x,E n )^+(x ,E n ) 

ra=0 



J dxm± t {x,xo) = I dxmX, t =o{x,x ) 



(A + 1) / dx 



1+A J dx'[tp(x',e,c) 



(20) 



(21) 



It means that the unique restriction to A is A ^ —1, that as it was early noticed corresponds 
to the absence of singularities at the quantum-mechanical level. It is hard to make an 
analogous study for mt t (x, xq) in a general case. However, for the Ornstein-Uhlenbeck 
process the normalization condition does not remove the parametric ambiguity at the 
definite choice of xq. 



3 The Ornstein-Uhlenbeck process 

Let us demonstrate the proposed scheme of getting new stochastic models with the well- 
known Ornstein-Uhlenbeck process. The Fokker-Planck equation which describes the 
Ornstein-Uhlenbeck leads to a quantum-mechanical potential with the harmonic oscillator 
potential. The factorization energy coincides with the ground state energy in the case 
with the harmonic oscillator potential. 

2 2 \ 

V- + ^- x A^e) = e V {x,e) (22) 

E n = nDuj,if)Z(x,E n ) = iu/D) 1 ^ H n [x^uj jD)t->\ n = 0, 1, ... (23) 
where e = (the energy counts from a). The potential entering the Langevin equation: 

UZ(x) = ux 2 /2 = D£ 2 /2, £ = ^jDx. (24) 

For the Ornstein-Uhlenbeck process it is easy to see that, at least for the case of 
x = 0, the normalized condition does also not fix the value of A. For an arbitrary 
stochastic process the same is hard to prove. Calculating m+ 4 (£,£ ) we use the wave 
functions ffTTl) . First of all it has to be pointed out of having the equilibrium value of 

<<-<«»> = < A + : » P \tvw ' m = 7?/ Ae " 2 (25) 

As it was mentioned in the above, the normalized condition of the probability density 
does not remove the A - parametric freedom. Using the wave functions ^/+(x, E n ), which 
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Figure 1: Dependence of m+ t (x, 0) on x 
and A at z = 1. 




Figure 3: Dependence of m% t {x, 0) on x 
and A at z = 0.25. 

are expressed in terms of *&Z(x, E n ), after 
expression for the probability density: 



Figure 2: Dependence of m\ t {x, 0) on x 
and A at z = 0. 




Figure 4: Dependence of m\ t {x, 0) on x 
and A at z = 0.5. 

simple but tedious calculations we get the 



m tt{£>£o) ^{nDj (l+A/2+A/2<J>(£)) 2 + 2 (d) d£ l+A/2+A/2S(£) X 

(1 + A/2 + A/2<&(*)) (l + $ (££)) - rf^' 2 [l + $ (^)] 

The form of 77i+t(C)Co) can be viewed from the following plots under the following 
parametrization ( = 0, z = e~ w *. 

Fig. ([Till]) show the results of numerical evaluation of the dependence of m+ t (x, 0) on 
x and A at different z. Thus, when z = 1 ml rt {x 1 0) has well pronounced 5-like shape in 
concordance with initial conditions. With decreasing of z m\. t (x, 0) changes its shape and 
this changes are especially significant at A — > — 1. In this case the shift of the maximum 
of probability density occurs, as well as emergence of substantial asymmetry. 
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4 Conclusions. 



The described procedure of the obtaining of exactly-solved stochastic models allows to use 
the results of numerous works concerning exactly-solved as well as quasi-exactly-solved 
quantum mechanical problems. The distinctive feature of this approach is an existence 
of parametric freedom in potentials, that enter the in the Langevin equation, as well as 
in transitional probability densities. This situation takes place at any modification in the 
spectrum of initial Hamiltonian. It's important to mention, that normalization condition 
for mi t (x,x ), most probably, could be fulfilled only when certain relation between x 
and A exists. This means that stochastic model with U*(x) is "bad". At the same time 
rri\_ t (x,xo) is a "good" transitional probability density, which preserve the parametric 
freedom. 

As is well known, [TJ E], local modifications of the shape of the potential in Langevin 
equation could lead to substantial changes of such quantities as times of the passing 
through barrier and time of live in metastable state. Although this quantities are mainly 
determined by first non-zero energy level in spectrum of H^, the A-freedom allows for 
significant variations of their values. This especially important for potentials with several 
local minima, which could emerge when constructing isospectral Hamiltonians with fac- 
torization energy e < E Q . In this case potentials with two wells emerge, symmetric, 
as well as asymmetric and A-freedom reveals in substantial modification of the shape and 
height of the barriers, what leads to changes in the rate of the inter-well transitions. 
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